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ABSTRACT 


The conversion of electromagnetic to acoustic energy is 
a broad and complicated subject. Some topics can be found 
in the literature but its thorough treatment is rather 
complex. 

Engineering approaches to the problem are however possi- 
ble through useful simplified models. Simplified treatments 
of the three major mechanisms leading to the production of 
mechanical forces on materials supporting an electromagnetic 
field are to be presented. These mechanisms are referred 
to as radiation pressure, transient surface heating and 
stimulated Brillouin scattering. 

A discussion of the experimental apparatus and procedures 
used in an attempt to experimentally verify some of the 


results given in the theoretical formulation, is also made. 
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LT. INTRODUCTION 


A) 6CU NATURE OF THE PROBLEM 

The generation of mechanical forces on materials 
supporting an electromagnetic field was first predicted 
several decades ago although its detailed theoretical 
treatment had been done much more recently. The literature 
concerning this phenomenon is not plentiful and is rather 
complex. Since the efficiency of the energy conversion is 
extremely low, the experimental verification of the theor- 
etical predictions was only possible after high power 
electromagnetic sources became available. Radiation from 
powerful radars as well as from high power pulsed lasers are 
the most suitable sources of electromagnetic energy for the 
study of the several mechanisms involved in this conversion 
process. 

The possible interpretations of the mechanical forces 
produced in a material medium subjected to interaction with 
an electromagnetic field, may be divided into three major 
types. That is, the forces may appear due to one or more 
of the following mechanisms: 

-(1)  yvradiation pressure, 
(1i) transient surface heating, 
(1ii) stimulated Brillouin scattering. 
It is a familiar fact that electromagnetic waves carry 


energy through the medium where they propagate in accordance 





with Poynting's vector. Less familiar, however, is the 
fact that electromagnetic waves may also transport linear 
momentum. 

Originally, the term radiation pressure was used to 
designate the pressure exerted on an object by shining 
light on it. This was predicted theoretically by Maxwell 
about 1870 but the first measurement of radiation pressure 
was made latter in 1901-1903. According to Maxwell's pre- 
dictions, if a parallel beam of light carrying the energy 
U falls on an object for a time t and is entirely absorbed, 
the magnitude of the momentum p delivered to it is given 
by p = U/c, where c is the speed of light. The direction 
of = is the direction of the incident beam. When the energy 
U is entirely reflected, the magnitude of the momentum de- 
livered is twice that given above, or p = 2U/c. On the other 
hand, if the energy U is partly reflected and partly absorbed, 
the delivered momentum will lie between U/c and 2U/c. 

The resultant forces are so small in relation to forces 
of the daily experience that they cannot be ordinarily 
noticed. 

Nowadays radiation pressure has a much broader signifi- 
cance. It is interpreted as the mechanism responsible for 
the production of forces in a medium due to a rate of change 
of electromagnetic momentum associated with moving energy 
as its equivalent mass. This rate of change of electro- 


magnetic momentum is the consequence of energy gradients 





established in the medium due to conversion of energy to 
heat or to progressive changes in the storage energy; it 
can also be the consequence of the existance of velocity 
gradients as it is the case when part of the incident 
energy is reflected back at the interface between two 
different media. 

The radiation pressure forces are in general much smaller 
than the forces produced through the transient heating of 
the surface of a medium. 

The generation of acoustic waves utilizing the thermal 
effects accompanying the absorption of electromagnetic 
radiation by an elastic medium was first experimentally 
verified in 1881 and since then just a few theoretical and 
experimental works have been developed. The subject seems 
to be receiving now more attention probably due to the fact 
that the recent advances in high power laser technology 
gives to this process the possibility of concrete and useful 
practical applications. 

When a pulse of electromagnetic radiation passes through 
an absorbing medium, its intensity is exponentially reduced 
producing a rise in the temperature of the medium. The 
thermal shock occuring as a consequence of the intense sur- 
face heating produces a thermal expansion originating an 
acoustic wave which propagates away from the heated surface. 
That is, a simple thermoacoustic source is established in 


the medium. 
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It is possible to extend this concept of a simple 
thermoacoustic source to the one of a thermoacoustic 
array which may be produced in the medium when it is sub- 
jected to a train of periodic pulses of high energy. [In 
this case, a highly directional acoustic wave is launched 
from the heated region and propagates in a direction normal 
to the electromagnetic beam. 

As suspected, the analysis of this thermodynamic con- 
version process involves the knowledge of theoretical prin- 
ciples related to fields like solid and fluid mechanics 
and thermodynamics. 

So far, the mechanisms leading to the conversion of 
electromagnetic to acoustic energy were described by an 
interaction process between an electromagnetic wave and the 
matter supporting it. Under certain conditions, sound may 
also result when two intersecting laser beams of very high 
energy interact with a medium. If the angle of interaction 
of the two beams is the appropriate, optically induced 
generation and amplification of sound will be verified. This 
is a particular case of the scattering phenomenon and is 
called stimulated Brillouin scattering. Note that this is 
not the same mechanism as the transient surface heating. 
While the latter utilizes the thermalization of modulated 
electromagnetic energy, the former is based on the mechanics 
of travelling wave parametric interactions where the resul- 
tant acoustic wave has a frequency equal to the difference 


in the frequencies of the two interacting light beams. 
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B. BACKGROUND CONCEPTS 

A more detailed theoretical analysis of the major 
mechanisms responsible for the conversion process will be 
presented next. But before that, it seems to be worthwhile 
to make a brief but hopefully clarifying comment about some 
unfamiliar terms belonging to fields like tensor calculus, 
mechanics and thermodynamics. The terms to be introduced 
here are mentioned in the theoretical formulation to be 
presented. 

Tensor analysis — this is a rather advanced subject in 
mathematics and so just a few basic concepts will be 
mentioned. 

In a three dimensional system consider two vectors A 
and B such that each rectangular component of vector B is 


a linear function of the components of vector A, that is 


to 
Il 
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If the association of the components of vector B with 
those of A is to be maintained for any rotation of the 
system coordinates, then the coefficients 7G have to trans- 
form in a specific manner. By definition, a tensor of rank 
two is a linear transformation of the components of a vector 


K into the components of a vector B which is invariant to 


HEZ 





rotations of the coordinate system. The nine components 
ace of the linear transformation are invariant scalars 
and are called the tensor components. 

The divergence of a tensor of second rank 1s a vector, 
or tensor of first rank, and PH divergence of a vector 
is an invariant scalar, or tensor of zero rank. The concept 
of a tensor can, of course, be generalized to include ten- 
sors of arbitrary rank. 

Stress~- suppose that a cut perpendicular to the x axis 
is made in a body in equilibrium. In order to maintain the 
equilibrium, an internal force would have to act in the 
exposed face. By definition, the stress at that point is 
the intensity of such a force, that is, the force per unit 
area. If AF is the resultant contribution of these internal 
forces on an element of area AA, then the components of AF 


along the three coordinate axis, AF Ly oc and AF are 


called stress components and are defined by 


| AF 
vx © os AA 
OF al 
xy ea MA ‘xz ee RA 
o., is the normal stress and ee co are the shear stresses 


and from its definitions it is clear that normal and shear 
stresses are, respectively, the intensity of force perpen- 


dicular and parallel to the plane of the area element. Normal 


eS 





stresses are termed tensile if they pull away from the cut; 
those pushing against the face are called compressive 
stresses. 

Obviously, any deformation of a body due to external 
forces will give rise to internal forces, or stresses. The 
stresses may then be thought of as entities describing the 
state of the internal action of a body. 

Imagine now, as an extension of this concept, that a 
cube of infinitesimal dimensions is isolated from a body 
and has its faces perpendicular to the coordinate axis. 

The stesses exposed by the cube are then the three normal 
stresses Oney oy, and oF and the six shear stresses T_, 


xY 


oe a i and t .. It can be shown [Ref. 1] that 
ny x 72 ZY ZX RZ 


T = T 7, T = T and T = Users 
XY yx eZ ZX a ZY 


The stress components may be arranged into a matrix 


format representing the state of stress at a point: 


O T ie 
x xy Ns 
T O T 
yx Y yz 
t v O 
ZX Zu Z, 


This is called a second rank stress tensor. 

Strain — the concept of strain is related to the 
deformations of a body. The normal strains in the axis 
directions, Ee Sy and eos are defined as the change in 


length per unit length of a line segment in the direction 
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under consideration. It follows from this definition that 
strain is a dimensionless quantity. Shear strains Ya5 
may be also defined; they represent changes in angles 

between line segments. The subscripts in Y53 have similar 
meanings to the subscripts for shear stresses. 

The strain components can also be assembled into a 
second rank strain tensor. 

It should be noted that stress and strain are closely 
related quantities. 

Elasticity — is the property of a material which enables 
1t to return to its original size and shape after the cause 
of a deformation has been removed. 

Bulk modulus of elasticity 'B' — also called modulus of 
dilatation, is defined as the ratio of the change in pressure 
on a body Ap, to the resulting fractional change in volume 
-AV/VV, or B = -VAp/AV. 

Specific heat 'C' — this is the heat capacity per unit 
mass of a body; heat capacity is the ratio of the heat AQ 
supplied to a body to its corresponding temperature rise 
iieesco, C = AO/mAT. 

Thermal conductivity 'K' — is a constant of proportion- 
ality relating the time rate of heat transfer dQ/dt across 
a certain area A and the corresponding temperature gradient 
at Gs worrgo/adt = -KA dT/dx. 

Mechanical equivalent of heat 'J' — since heat is just 
another form of energy, any energy unit could be a heat 


unit. Hence, the mechanical energy equivalent of heat 


5 





energy, that is, the number of joules equivalent to one 
calorie is called mechanical equivalent of heat, 
J = 4.185 J/fcal. 


Coefficient of linear thermal expansion 'a' — the 
change in any linear dimension of a body is called a linear 
expansion. This coefficient is then a constant of propor- 
tionality relating the temperature change AT of a body, 
its original linear dimension 2 and the corresponding change 
imme his dimension AL, or AL = alAT. 

Adiabatic process — is a process taking place in such 
a way that no heat flows into or out of the system. In 
practice, since the flow of heat is somewhat slow, any 


process can be made approximately adiabatic if it is per- 


formed quickly enough. 
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Mee LHEORY 


A. RADIATION PRESSURE 
1. Elastic Stress Tensor 

The study and description of a certain category of 
physical phenomena like for example, volume deformations, 
camnot be done correctly by the use of the scalar or vector 
concepts. The adequate treatment of this type of problems 
requires the use of the tensor concept; tensor calculus, 
a somewhat unfamiliar tool, plays then an important role 
in the analysis of such problems. 

When dealing with mechanical forces on a material 
medium, it is necessary to introduce an entity called elas- 
tic stress tensor. This is a tensor of second rank and so 
it has nine components, the stresses, representing the 
forces exerted on unit elements of area. It turns out that 
the divergence of this tensor is a vector, or tensor of 
first rank, giving the force per unit volume caused by the 
elastic stresses. 

JeleecErOnmagneele Stress Tensor 

An electromagnetic field is now considered; as 
usual in electromagnetic theory, its definition is based 
upon the four vectors E, B, D and i satisfying Maxwell's 
equations. E and H are, respectively, the electric and 
magnetic field vectors. The electric and magnetic flux 


density vectors, D and i are related to EF and H by the 


yl 





-> 
constitutive equations D = cE and B = WH where € is the 


permittivity and uw is the permeability of the medium under 
consideration. 

The existence of an electromagnetic field ina 
material medium leads to the definition of the so called 


electromagnetic stress tensor; its nine components Ty 


x Txy! 
ee! Tx! etc. The electromagnetic stresses are defined 
[Ref. 2] by the following relations 

i ie 2 cee 2 ee 
ie SE E EO) +5 (Hy oe He) 
Similarly for Ty! ee 
1 = T = Cr F +) wn H 
xY pes x Y x Y 
Sais eines Otnash 2! ce (ae), 


The divergence of this tensor, div T, is a vector 


quantity whose components in carteSian coordinates are 


obtained according to the rule 


oT ey oT oT Ls 
eS a ty se 


Simibariny £or (div TY and (div ae 


(2) 
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Using Maxwell's equations it can be shown [Ref. 3] 
that 

div T = oF + oxB + ese (ExB) (3) 
where o and J are, respectively, volume charge density and 
current density. 

As already mentioned, the divergence of the elastic 
stress tensor is interpreted as a force per unit volume; this 
important result can be extended to the case of the electro- 
magnetic stress tensor so that its divergence may be inter- 
preted as an alectromagnetic force per unit volume. 

3. Electromagnetic Momentum 

It is apparent that the terms oF and xB in the 
right-hand side of eq. (3) have dimensions of a force per 
unit volume. For a medium with o and ci Pai OUG meu. cd 
volume V, these quantities represent then the forces exerted 
on unit volume elements of charge and current, that is, they 
are the well known electric and magnetic forces exerted on 
charges and currents. This means that a portion of the total 
force acting on the unit volume appears as a force on the 
matter within that volume. Therefore, the remaining force 
exerted on the unit volume should be associated with the 


terme d/dt (ExB) ; 


ee 


As the dimensions of B = WH are MQ -T ~ and those 


Z 


of D = cE are QL , the quantity  (ExB) is dimensionally 


eo 
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il 


on *Mo7 T -~ = ML “T ~ = MLT 


volume since the dimensions of momentun, D = mv, are mur 2. 
As a consequence, the term e (EXB) 1s interpreted as the 
electromagnetic momentum density. Now it is clear that the 
remaining force exerted on the unit volume by the external 
field acts by increasing the electromagnetic momentum within 
that volume. 

The precise manner in which electromagnetic momentum 
density should be defined has found some differences of 
Opinion. The existence of a duality of treatment of such 
concept results mainly from the two different ways in which 
matter can be visualized. On the one hand, a macroscopic 
presentation of a medium may be adopted when dealing with 
electromagnetic propagation in that medium; on the other, 

a microscopic picture describing matter as an arrangement of 
dipoles, doublets and charges surrounded by vacuum, exists 
associated with the concept of energy whose value depends on 
the relative velocity, but which always moves, regardless 
of the medium, with the velocity of light in vacuum. 

According to the latter interpretation, the electro- 
magnetic momentum density is defined [Ref. 4] by Moe, (EXE) , 
where c = Re is the velocity of light in free space. 

To summarize, it can be stated that associated with 
an electromagnetic field, there is a momentum distributed 
with a density Dp which is given, according to the specific 


interpretation, by one of the two possible expressions 
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Ue, (ExH) (4) 


ue (EXE) (5) 


Or 
ll 


Throughout this work the latter interpretation will be 
mostly used. 
syeeioeces inthe Electromagnetic Field 
The divergence theorem, a useful mathematical 


relation in electromagnetic theory, 


6 D-aS = V°D dv (6) 
S 


f 
V 
states that the integral of the normal component of any 
vector field over a closed surface S, is equal to the 
integral of the divergence of this vector field throughout 
the volume enclosed by the closed surface. 
This vector divergence theorem has a tensor analogue 
Stating that the integral of the divergence of a tensor 
throughout a volume V is equal to the integral of its 
normal component over the surface S enclosed by V. 
Integrating eq. (3) over the volume V and applying 


this theorem, the result would be 


> + > > 0 > > 
T-dS = f pE AV + f (JxB)dv ea { (EXB) dv 
V 


V V 7) 
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Consider for the moment a stationary distribution 
of charge and current; the fields are independent of time 
and so the last term on the right of eq. (7) is zero. The 


other two terms 


[SpE aV=F, , f (xB) dv = FL (8) 
represent the net mechanical forces acting on the volume 
distributions of charge and current, where o and J are the 
charge and current densities within the closed surface S. 

The total force transmitted by the electromagnetic 
field across the surface S or, the resultant force exerted 


on the charged matter within S is 


Fo+F = f (o£ + JxB) av (9) 
V 


It can be shown [Ref. 3] that 


Fo= sf pEav = f f[e(E-h)E - § £°n) as (10) 
Vv S 
a ~ oe i fee oe 
BF, = f (xB) dv = sf f{ = & n)B s BYn] ds (11) 
S 


where n is the outward unit normal at a point in the surface 


enclosing V, EB? = ES *F EY ' ES = jz | and 
eases = iB \- . 
x 4 Z 
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By Newton's second law 
d 
auaaet a= so P (12) 


os is the total linear momentum of the charged matter within S. 
Consider now that the volume V contains no charges 

Or currents. If the field is stationary, the net force 

across S is zero. But if the field is variable, the third 

term in the right-hand side of eq. (7) is not zero and so 


this equation becomes 


Kaecicue= = f e(BxB) av (13) 
S Vv 


The electromagnetic momentum concept shows up 
now to help in understanding this apparent action of a force 
arising in neutral matter immersed in an electromagnetic 
field. As stated before, c (ExB) = D 1s the density of the 
electromagnetic momentum associated with the electromagnetic 


field. The integral 


Sf pav = B, (14) 
V 
1s the total momentum of the field contained within V and 
then, as eq. (13) says, the force transmitted across S goes 
to increase the electromagnetic momentum of the field 


Cinthia Vis 


a 





Consider at last a system composed of charges and 
variable fields within a bounded region. The conservation 
of momentum theorem is, therefore, according to eq. (7) 


expressed by 
f. (@ +B) = 7s tad (15) 
iis m © 
hwecollows trom eqs. (9), (10) and (11) that 


(eE* + vH*) a] as (16) 


bof ke 


Fae eee Gel (i) acto (Hn) He 
S 


This is the resultant force on the charge, current and matter 
within the surface S. 

Equation (16) is a valid relation for both stationary 
and dynamic situations and the validity of its application 
to the case of variable fields can be justified if basic 
Principles of the theory of relativity are invoked. Then, 
the rate of change of both mechanical momentum of the matter 
within S and electromagnetic momentum of the field within V, 
may be equated to the right-hand side of eq. (16) which must 
now be interpreted as the inward flow of momentum through 


the surface S per unit of time [Ref. 3] 


8 ye) ds 
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Transforming this surface integral into a volume 
integral and using Maxwell's equations, the resultant 
expression would have the form 


Z 


( Ye +H Vu) +22(5 xB) Jav 


(18) 


Nol 


(Pp +P) = S[pEt+taxB- 
Vv 


iiomenenecsloneq@vesmithie net LOotal force acting on 
the volume V; it 1s equal to the rate of change of the 
total momentum, electromagnetic and mechanical. 

The forces exerted by the field on charges and 
neutral matter give rise to an increase in the mechanical 
momentum so that F = d/dt P This force can be determined 
1£f the right-hand side of Eq. (18) can be split into two 
parts to be identified with E and Poe However, the analysis 
of the problem may become difficult since this separation 
is not an obvious task to carry out. That is why the 
literature on the subject shows that several hypotheses have 
been suggested to define what part of the total momentum 
1s associated with the strictive forces in the medium pro- 
ducing mechanical deformation. 

In order to simplify the analysis of the mechanical 
forces arising from the interaction of a travelling electo- 
magnetic wave with a material medium, it will be assumed 
that the medium is homogeneous, devoid of free-charges 
and of mechanical elasticity, and that steady-state condi- 


tions apply. 
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Under these assumptions, Eq. (18) reduces to an 
expression for the rate of change of electromagnetic 
momentum giving two force components which are represented 
by the terms 3 xB and 3/at(D xB). As 5 xB is independent 
of time variations, it may be interpreted as the force 
that would arise in an equivalent DC field; on the other 
hand, 9/3t(D xB) can be interpreted as the force associated 
with the travelling electromagnetic wave. 

It is then clear that to calculate the mechanical 
forces on the material medium, the time rate of change of 


the electromagnetic momentum has to be known. 


Bee fouces Duetto Electromagnetic Energy and Velocity 


Gradients in a Material Medium 
The total energy density due to the electric and 


7 
magnetic fields is given by 
W= (eE* + uH*) (19) 


This expression represents the instantaneous energy stored 
in the electromagnetic field per unit area. According to 
the theory of relativity, mass m and energy W are equiva- 


lent and must obey the relation 
_ Ww 
SS gy 2.09) 


Also, if mass m moving with velocity v has a rest mass 


M+ the following relation must be valid: 
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m 5 
m= pe (21) 


/1 = (v/e)* 


The binomial expansion of the right-hand side of Eq. (21) 
will result in the expression 


: - a ea 
a a: - mM a mls a! oe: aa rk (22) 


ms the kinetic mass, is a function of the velocity v; 
it vanishes when v = 0. 

Assuming that Eq. (20) can be applied to moving 
electromagnetic energy and its equivalent mass, then the 
energy W can be split into a rest energy W, and an energy 
of relative motion W 


k 


W = ae = (m_+m,)c (23) 


Consider now a uniform plane wave with an angular 
frequency w and a phase shift constant 8 in a medium with 
conductivity o. The phase velocity v of the travelling 
wave is given by 


1/2 


ices Ieeeintey ie + 1 


If the wave comprises an energy density Tei = (cE + uH*) , 


then the rate of flow of energy per unit area is: 
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Sue = eae jy = Wiv = EH (25) 


The magnitude of the electromagnetic momentum 
lel > > : ; 
density p = DxB = we(ExH) associated with the moving 


energy W. 1s 


W_v 

= = EH — BH — a 
p = WEEH = w/e a pee C216) 

O O 


where Me ley el 
The mass-energy equivalence principle applied to 
the energy We would be expressed by the relation m= Wo/ve. 
Note that the energies ns and Wy must not be considered 
the same quantities. Whe 1S an additional energy of motion 
for a rest energy we and it translates to equivalent mass 
in terms of . 
For a transverse ee eeenaa netic wave progressing 
in the +Z direction, the electromagnetic momentum associated 


Wlth an element of thickness 6Z and unit cross-sectional 


area is given by the expression 


UU 
V 
O 


As force is the time rate of change of momentum, the 


corresponding elementary force would be 


ie Sy eee aaa QV) 5¢ (28) 
dt + oP Se ce a dt 


O 
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From basic calculus 


so that 


2 OW WV 


sF= (Uy —3 + 3 $2 (29) 
Mie OZ Mies 


The integration of Eq. (29) allows then the computation 

of the total force per unit surface area and if the time 
averages of the force, F, and energy, W, are considered, 
the expression for the total time average force per unit 


surface area will be 


OW 


2 
a a L os dV 
F = U2 ) az + 2 {Wav a az | (30) 
O O 


Eq. (30) denotes the presence of energy and velocity 
gradients and since mechanical forces can only operate on a 
material medium, these energy gradients must arise either 
from dissipation producing heat or from a progressive 
Space change in the ooeeed energy. Also, a space change 
in the propagation velocity yields a force component as it 
is the case when power is reflected at an interface between 
two different media. 

As stated before, the gradients of energy in a 


source-free medium resulting in a rate of change of the 


OS, 





electromagnetic momentum, may be produced in two distinct 
ways: 

(i) by conversion of energy to heat, 

(11) by transfer of stored energy to or from the 

material medium. 

In the source-free conditions assumed, any additional 
energy put into the material as the wave propagates through 
it, whether to provide for heat losses or a change in the 
stored energy, must come from the empty space which the 
material occupies since conservation of energy requires a 
constant value of its total content. Any such transfer of 
energy from space to material medium is associated with a 
force acting in the +Z direction, that is, 1t 1S a positive 
force. Then, the dissipation of energy as heat or a 
progressive increase of the stored energy must produce a 
positive force. | 

The resultant force produced in a material medium 
through the process described is called radiation pressure. 

Palen eoradlents)— A More Detarled Analysis 
For the conditions assumed, it follows from 


Maxwell's equations that 


gE JB 
Ze sot 3) 
dH 3D. 
= 7 =e E.. =f a MSZ) 
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Also, the expression a = (cE~ + LHS) may be written 


as 


Saal 2 2 
Ne = ZleEL a uA.) (63a) 
Combining Eqs. (31) and (32) with the partial derivative 
of Eq. (33) with respect to Z, an expression for the 


gradient of the total energy in the electromagnetic field 


can be readily obtained 


aw. 


ome 2 ee 
Van ag Bee a Qui) 
dB aD 
= - = yb —X lope 9€ , 4? Ou 
Hoe A cE, se uH ce x (ES 5 a 1) 
(34) 


Under the assumptions being considered here, i.e., charge- 
free medium and unit cross-sectional area, it should be 
Poveda that the right-hand side of Eq. (34) corresponds to 
Eq. (18) and therefore it gives the rate of change of the 
total momentum including the part arising from strictive 
forces producing mechanical deformation which, however, is 
assumed to be negligible in the present discussion. 

Since the medium is homogeneous the last term on 
the right side of Eq. (34) vanishes. 

When the fields have a sinusoidal time variation, 
which is normally the case, the time average of the energy 


gradient corresponding to Eq. (34) is easily computed using 


Si 





complex wave functions. Thus, representing the field 
components as complex amplitudes including both the phase 
factor and the space dependence of the wave function, the 


result is 


JE 


dH 
bed - *) — ey * x 
= HONE Coa weRe (EY se ) peRe (He se ) (35) 








where Re and the asterisk stand, respectively, for the 
real part and complex conjugate of the quantities. But 
noting that 


aH, JE, 
~ aera * 
Re (EY ) ee ) (36) 








ac dt 


it happens that Eq. (35) reduces to the form 





= -yoRe (5 E. HY) 37} 


This time average of the gradient of the total energy per 
unit volume is associated with conversion to heat, a fact 
that is a consequence of the finite conductivity o of the 
medium; as this gradient is a negative quantity, it repre- 
sents a progressive transfer of energy to the material with 
a gradient +SuoRe (E,, 7 and so the corresponding force is 
positive. 


Eq. (37) does not reveal any gradient of stored 


energy in the material; therefore it is necessary for 


a2 





further analysis to show the possibility of the existence 
of such a gradient in the medium responsible for a 
supplementary force. 

Representing the permittivity and permeability of 
free space respectively by EG and Woe the instantaneous 
energy stored in the material medium alone per unit 


volume 1S given by 
ee Z a 2 
Wi2 = 5fe ES) EY + Cy oo (38) 


Differentiating Eq. (38) with respect to Z, taking its 
time average uSing complex notation and assuming a homo- 
geneous medium with u = Woe the final expression for the 
gradient of stored energy would be 


JE 


= ule -e,)Re(s He x) (39) 
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In the case of a pure travelling wave this gradient will 
be negative [Ref. 5] because as the field progresses 
through the medium it attenuates due to its finite conduc- 
tivity. The gradient represents then energy progressively 
coming out of the medium and the force produced will be 
negative. 

So, in general, it appears that a force component 
may result from a rate of change of electromagnetic 


momentum associated with stored energy. 
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7. Calculation of Radiation Pressure in an Infinite 


Homogeneous Medium 


Eq. (30) is the basic equation for the calculation 
of the forces produced in a medium with characteristics 
u, € and o, due to a rate of change of electromagnetic 
momentum associated with a travelling wave moving in the 
positive Z direction with a constant velocity v. Since 
the velocity is constant, the second term on the right- 
hand side of Eq. (30) is zero and so this equation becomes 


for a thickness d, 
Z 
= V 
Vv 
oO 


aw _/a2 can be decomposed into two gradients of energy 

ie 
according to Eqs. (37) and (38) so that, using the appro- 
priate signs to yield positive forces, aW_/32 can be written 


as 


aw, oW oW 








ae al a2 
Piz az CO ene 
and Eq. (40) becomes 
FFL t+ Fs, 
2 da dE 
av i * = ee x 

a 2 pee ee E_) Re (FEY ar ) ] dz (42) 
O 
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where Fy and ie are the average forces resulting from the 
gradients of energy associated, respectively, with dissi- 
pation as heat and with stored energy. 
a. Forces in a Medium Extending from Z=0 to Z=~ 
For a finite conductivity medium extending from 
Z=0 to Z=~e, the total power associated with a wave is 
absorbed which means that the overall stored energy is 


zero. Thus, no net component of force F. can be expected. 


2 


When this situation occurs, Eq. (42) takes the form 


= =m _v = a a 


oO 


The quantity Re (= E, me) is the time average value of the 
instantaneous power density s carried by the wave which 
attenuates with Z at a rate determined by a, the attenuation 
coefficient of the medium; it can be shown [Ref. 6] that 


Ss is given by 


Z 
EB 
S Re (= Ey 7 Pare os § e sé (44) 
where E. is the value of E, at Z = 0 and In| is the magni- 


tude of the complex intrinsic impedance of the medium 


n= |n| fe. 


Then the force coming from energy dissipation 


is given by 


2 co 2 Ss 
= = _v — -20Z _V x 
a Fy teeny 7 sf wos,e az —xHO (55) (45) 
Vo 0 O 
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and since v = 2a/u0 


F, = ae: S$) (46) 
O 
b. Forces Due to Reflection in the Medium 
If the wave being propagated through the medium 

finds a partially reflecting surface, then a certain per- 
centage s = Wv of the power it carries and its associated 
fuecerOmagnetic momentum, suffer a complete reversal and 
so the second term on the right side of Eq. (30) is no 
longer zero. This means that a force of reflection will 


be produced 


+V 


~- 4 — 9V _2sv 
emma a 57 OF se 
Yo aah’ s Yo 


For the particular case of a lossless medium (o = 0), 
the fact that v = a allows a further simplification of 


Eq. (47) and since y = eee and «=e ¢€ 





O Te 
= — We 8 25 
"3 ee YEy Hy Oe 


c is the velocity of light in free space, and i are 
respectively the dielectric constant and relative permea- 
bility of the medium outside the reflecting surface. 

As seen from Eq. (48), the force due to reflec- 
tion is directly proportional to ve. , the index of refrac- 


tion of the medium. 
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8. Rate of Change of Electromagnetic Momentum for 
A DC Field and its Corresponding Force 


When frequency approaches zero, the velocity v 
approaches zero too and the same happens to the electro- 
magnetic momentum as can be seen from Eq. (26), for example. 

Consider the material medium immersed in a DC 
field. If a rest energy Wy analogous to W- is considered, 
it may be translated to equivalent mass ™m, using the 
MeLOCity c. 

As a logical extension to this case of 'DC flow 


of energy', momentum may be then expressed by the relation 


he 


"bp 
Bea ee = yey 7 (49) 
So, even in the case of a DC field, the concept of an 
electromagnetic momentum responsible, when absorbed by the 
medium, for a mechanical force, is still valid. Without 
gOing in much detail, it seems natural to accept the idea 
of a smooth transition from wave propagation to DC trans- 
mission. It should be noted that this electromagnetic 
momentum density concept is equivalent to the one described 
Byena.. (4). 

It will be shown next that under DC conditions, the 
production of force on the material medium via conversion 
of electrical energy to heat is still maintained. 

Wy/¢ represents the volume density of electromagnetic 


momentum for a DC field in a homogeneous medium. When the 
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field components are transverse to the direction Z, the 
electromagnetic momentum associated with an element of 


thickness 6Z and unit cross-sectional area is 
W 
Sp = = $2 (50) 


so that the corresponding elementary force is 


W dw 
cd _ b ~ b 
By integration, the total force per unit surface area is 


given by 


For time-invariant field components, the displace- 
ment current density 9D /at vanishes and so Eq. (32) 


becomes 


0H 
VELL G EY C53) 


Using Eqs. (31) and (53), Eq. (34) reduces to 


OW, 
77 a Bose (54) 
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For a thickness d of the medium, the force per 
unit area arising from dissipation of power, according to 
Bae) (52), 1S given by 

d 
Eb Sato EH, dz (55) 
0 

Obviously, no gradient of stored energy contri- 

butes to this mechanism of production of force. 
9. Summary of Important Results 

The discussion of the forces arising in a material 
medium under the action of an electromagnetic field, was 
carried out in a simplified version and results ina few 
important relations. 

It was shown that the process responsible for the 
production of such mechanical forces, is a rate of change 
of electromagnetic momentum associated with moving energy 
Or its equivalent mass. 

Three different mechanisms may be concurrent in the 
production of the resultant force, each one being, however, 
a consequence of distinct causes originating a change in 
momentum. These three major mechanisms described before 
are repeated here along with the corresponding equations. 

a. Change in Momentum Due to Energy Gradients 

Energy gradients in the medium may be produced 
in two distinct ways giving rise to different force 


components: 
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(1) By conversion of energa co heat 


= v* d i 
= ee) * 
O 


a) EVmenansiter Ofgstored energy to or from 


the medium 


= v4 d i OB. 
= = eae 
F ze “od u(e SO) a Dae xe ) daZ (57) 
Oo 


b. Change in Momentum Due to Velocity Gradients 
A force component can also arise if, fora 
given power (Wov), there 1s a change in propagation 


velocity 


mn al = 
- = a, ae 
V 
O 


az (58) 
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B. TRANSIENT SURFACE HEATING 
People Thermoacoustic Source 

Another process of the conversion of electromagnetic 
to acoustic energy to be described now, is the so-called 
thermal stress generation by transient surface heating of 
an elastic medium. 

This is a known process since the end of the 19th 
century but only recently, due to improvements in high power 
laser technology, it was possible to make an experimental 


verification of the complex theory governing this mechanism. 
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The solution of the related problems require a 
certain number of assumptions so that the mathematical 
models used apply only to the specified conditions. 

The transient heating of the surface of a body may 
be produced in several ways, namely, by electron bombard- 
ment or by electromagnetic energy absorption. The latter 
process is the one. to be considered here, that is, the 
conversion to an elastic wave of that portion of the 
incident electromagnetic energy absorbed by the medium. 

It will be adopted a one-dimensional treatment 
of the problem and assumed an uniform heating at and near 
the surface of a semi-infinite absorbing medium as well 
aS an incident electromagnetic intensity uniformly distributed. 

Under these conditions, a temperature gradient 
normal to the surface will be produced. This gradient is 
characterized by the physical properties of the medium and 
by the nature of the incident flux. The consequent ther- 
mal expansion gives rise to a production of strains in the 
medium which, in turn, leads tothe generation of stress 
waves propagating away from the heated surface. 

The time profile of the stresses and the amplitude 
of the waves produced for a given absorbed power density, 
depend upon the acoustic boundary conditions at the illum- 
inated interface, that is, upon the elastic constraints 
applied to the heated surface. These boundary conditions 


may, in general, be divided into two types: stress-free 
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and constrained by contact with another body. The stresses 
and wave amplitudes obtained in the stress-free condition 
may be smaller than those obtained under constrained 
heated surface. 

2. Equation of Motion 

Consider that the plane x = 0 in a cartesian 

coordinate system is the surface of the semi-infinite 
elastic medium being illuminated with uniformly distributed 
electromagnetic energy. Part of this energy will be re- 
flected at the interface but the remaining will penetrate 
the medium. If a and Io represent, respectively, the 
attenuation constant of the medium and the radiation inten- 
sity at the surface, then the radiation intensity I(x) at 
a distance x from the plane x = 0 is given by the familiar 


exponential law 


I(x) =I. ee (59) 


This means that the electromagnetic wave is being absorbed 
by the medium as it progresses through it. 

This absorption of energy causes a temperature rise 
at and near the surface and, as a result of this effect, 


a strain 25 ls produced 


Pou, c) 


= a = aO(x,t) (60) 
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where u(x,t) 1s the x component of particle displacement, 
ais the coefficient of linear thermal expansion and 
@(x,t) is the temperature rise above the uniform initial 


temperature. Recalling that the problem is being treated 


Taeone dimension 


CMe) me oulx,t) .. 
<a eon ae 


therefore the strains ST and Ey are also zero. 
In order to produce the strain EL in the absence 


of heating, a stress o.. would have been necessary 
o = — Bad (62) 


where B is the bulk modulus of elasticity. 

In the presence of both heating and stress, it 
can be shown [Ref. 7] that the corresponding stress-strain 
relationship for the case of negligible lateral inertia 


and shear, exhibited by liquid media for example, is given 


by 
¢ = Be, Baé B a Baé (63) 


The resulting equation of motion for the semi- 


infinite body is then expressed in the form 
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0 aa = igs 
yee ox 
See ou _ 
ae (B ax Baé) (64) 


where 0 is the density of the medium. Neotunge that the 
modulus of elasticity and the velocity v of the compressional 


wave propagation in the medium, are related by 


B = pv (65) 
eq. (64) becomes 
1 ou _ 3-u a eG. (66) 
ye te 337 22 


This is the nonhomogeneous wave equation to be solved for 
the elastic particle displacement. 
See Lenperacure Distributions 

In order to solve eq. (66), the temperature distri- 
bution function §(x,t) has to be known. 

For the case of an electromagnetic wave, if 1t was 
assumed that the whole energy was absorbed in the plane 
x = 0, relatively simple approximate temperature distributions 
could be obtained. However, absorption of energy actually 
takes place in a region of finite thickness which is related 
to the skin depth. If this fact is taken into account, the 
temperature distributions are obtained through a more 


difficult analysis but they will be accurate. 
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A particular temperature distribution is examined 
here: it is the one resulting from a constant input electro- 
magnetic intensity of short time duration, suffering an 
exponential attenuation in the medium. It is assumed that 
most of the input energy degenerates to heat and that the 
process is adiabatic. 

If radiation of intensity I, hits the interface in 
the time range t = 0 to t = T, a heat source in the region 
x > 0 will be created having a heat production rate deter- 
mined by eq. (59). The resultant heat flow can be considered 


1 1s much 


one-dimensional as long as the penetration depth a 
less than the diameter d of the surface receiving heat. 

A solution for this temperature distribution, where 
the heat flow is supposed to be in the positive x direction, 
was achieved by Carslaw and Jaeger [Ref. 8] but the final form 
of the function §6(x,t) is a rather complicated expression. 

Under certain conditions, the solution just mentioned 
may be simplified. In fact, an approximation can be found 
[Ref. 7] for the case of intense heating of media having low 
thermal conductivity, such as liquids, and for durations of 


the radiation pulses from 10 to 50 nsec. This approximation 


1s expressed by 


—-QAX 
at te 


7 =e oe 


Osa, it) 


where J and C are, respectively, the mechanical equivalent 


of heat and the specific heat of the medium. 
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This type of heating can be produced using, for 
example, a Q-switched laser with output in the range 
50-100 MW/om*. The temperature rises so obtained, although 
rapid, do not exceed tens of degrees but the experienced 
temperature gradients would be extremely large. As an 
example [Ref. 7], it could be expected in water colored by 


Zi 1 


a dye to give a = 10“cm ~, temperature gradients as high 


: °Cc/com and temperature rates of change of the order 


biome) «x10 
of 2.5x10° °C/sec at the surface during the laser pulse. 

If nonmetalic media is considered, the temperature 
decay for t > T is a slowly varying function of time, 
becoming Significant only for t greater than milliseconds. 
Also, the elastic tranSients produced through this process, 


propagate away from the heated place in periods of time as 


short as microseconds. So, it can be assumed that for t > T 


a, we 
Q@(x,t) = a Cue (68) 
Equations (67) and (68) represent then the source 
functions to be used in eq. (66) in order to determine the 


equations for displacement and stress. 


4, Solution of the Equation of Motion 


For Two Boundary Conditions 


a. Constrained surface (rigid boundary) 
Substituting eqs. (67) and (68) in eq. (66), 
the resultant wave equations during and after the heating 


pulse are 
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5 > vWte ; eect; 
2 
a = V a3 ae (69) 
at ox y-wre o% tog 
where 
ad I, 
W = ccm (70) 
The solution of these equations for the initial 
conditions 


23 = 11 GeO) = 0 (71) 


and the appropriate boundary condition for a constrained 


surface 
wo, tc) = 0 (72) 


would give the expressions for the particle displacement 
Weweeiwerie= time ranges T < £ < x/v , xX/v < t < x/v + T 
and t > x/v + T [Ref. 7]. Once u(x,t) is known, du/dx 
can be substituted in eq. (63) to obtain the expression 
for the stress in the medium as a function of time and 
distance from the surface. Since the condition ee <oerd 
1s met, an elastic plane wave is then launched from the 


surface and propagates away from it. 
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Again, due to the extreme complexity of the equations 
for displacement and stress actually obtained, it is desirable 
to try some sort of valid simplification to render them more 
Suitable for specific applications. This is possible if the 
values of x are chosen such that ax >>1, that is, the points 
x Should be those far from the source. Then, introducing 
a reduced time variable T = t - x/v, the approximate 


expressions for the stress would take the form 


—enn7 1 OL wicl 





(l-e je re ae 
val = = 
o(x,t) = => . [2 -e ONS a a O<T<T, (73) 
ele = ive ae ee 


From these equations it can be observed that the 


maximum stress occurs at T = T/2 and its value is given by 


val 
4 oO avi 2 
Onax ~ BIG al € ) (74) 
which approaches the limiting value of val ,/2dC BOR (yee oo les 
b. Free surface (pressure-release boundary) 


The boundary condition at the stress-free surface 


is now 


eum 0 ,t) 


ee = Bage0,t) = 0 (7) 


So, ee — 8 


or, equivalently, 
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omit et) = ag@(0,t) = | (76) 
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Under this condition and following a procedure similar to 
the one used for the constrained surface, it is possible to 
get a set of ee omianed equations describing the stress at 
observation points far from the heated surface, ax >> l, 


and in terms ofT 





(ues ae oe aaa 
val 
o(x,t) = an : [e a ene a , O<TST , 7) 
d= Be ae rT 


These equations tell that maximum stresses having the value 


val 
eS ere OL 
a cc ee ae 
appear at tT = 0 and Tt = T and approach a limiting value of 


val /2uc BOerne os 1, . 
5. Stress Waves in the Medium 
Solving eqs. (73) and (77) for various values of 
the parameter avT, a plot of stress versus reduced time can 


be obtained. 
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Typical curves obtained in this manner would look 


like those shown in Fig. l. 


‘stress d 





(a) (b) 


Fig. 1. Typical stress-time plots for constrained 
surface (a) and stress-free surface (b) 


The constrained surface situation gives rise to stress 
consisting of a unidirectional compressive pulse, while in 
the free surface case the stress suffers an abrupt change 
from compression to tension during the time T. 

To give an idea about the theoretical magnitude of 
the resultant stresses, the use of the free-surface expressions 
for the case of a 100 MW, 30 nsec laser pulse, illuminating 
water with a = 100 en would lead to produced stresses with 
peaks of +20 atm. 

As expected from the theory, the peak amplitude of 
the pulse decreases and its duration increases as the 


attenuation constant of the medium is reduced. 
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High power laser pulses may be then responsible 
for the production of thermal effects in media. The rapid 
changes in temperature occuring at a surface absorbing 
rapidly changing amounts of energy is due to the fact that 
thermal conduction is a slow process. The elastic waves 
detected at points far from the surface of a medium as a 
result of the absorption of an electromagnetic pulse of ~* 
high frequency, will have a frequency content consisting 
primarily of frequency components associated with the pulse 
envelope rather than with the carrier. 

A few experiments on the generation of acoustic 
Signals in liquids by laser pulses, had been performed 
recently and the results seem to confirm the theory behind 
the transient heating process of elastic wave generation. 

E.F. Carome, N.A. Clark and C.E. Moeller [Ref. 9], 
uSing a laser beam from a 0.1 Joules, 50 nsec single pulsed 
ruby laser, verified the validity of eqs. (73) and (77) for 
the case of transient heating of water surface. The measure- 
ments were performed with various concentrations of Prussian 
blue in order to change the attenuation constant of the water. 
The shape of the pulses obtained, closely approximate those 
predicted by theory. As observed, an increase in a, that is, 
an increase in the absorptivity of the medium, increased the 
elastic wave amplitude. This result indicates that the 
produced elastic waves could not have been generated through 
a process of radiation pressure by reflection at the surface 
Since an increase in absorptivity decreases the amount of 


reflected energy. 


ay 


6. Bienes Ney of Energy Conversion 


The conversion efficiency N of this process is 


defined as the ratio of the energy of the propagated elastic 


wave ES to the energy of the incident electromagnetic wave E 


Z 

ll 
‘alo 
m | 


(79) 


As may be guessed, this process is characterized by 


a very low efficiency. 


For both free and constrained surface conditions, the 


energy ES is a function of the area through which the elastic 


wave is propagated and the area over which the electromagnetic 


wave is incident. Assuming that both areas have the same 


value, the resultant expressions for N [Ref. 


fo 
case of free-surface 


where 


e avt aval 


ee - avTe 


and for the case of constrained surface 


va*I 
N = ——s- Fo (avT) 
XE oi 2 


a2 


for the 
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where 


= Oh fels —ovy l 
a ove + 3e ee! VL ea 3 


In these expressions oe 1s in watts per square centimeter 
while all’theother quantities are in the cgs system of units. 
A plot of the functions Fy and FE. versus avT, allow 
an easier and quicker numerical calculation of N. These 
plots show that the maximum value of Fy fo20. 5 foneavyt = 2.0 


and that F, approaches a maximum of 2.0 for large values of 


2 


CE. « 
As an illustration of a numerical calculation of the 


efficiency showing its dependence onthe properaties of the 


medium, the value N = eee se 


Z 


IF, (av) 1s obtained for 
water and N = 6.14x10 TGF, (ov) LOtecasbon tetrachloride, 
both for the free surface condition. 

A final remark about the functional dependence of N 
on the variable Ty should be made. From the anaysis of 
eqs. (80) and (82) it seems apparent that efficiency can be 
increased without limit as Ly 1s increased. However, this 
cannot be a true statement since beyond a certain value of 
input power the assumed model of stress generation is no 
longer valid. Some limitations of this model result from 
the following facts. First, the thermal characteristics of 
the medium impose limits on the radiation intensity since 
phenomena like boiling in liquids or melting and vaporization 


in solids may occur. Stresses still continue to be generated 
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although by a different process. Secondly, since at high 
intensities the thermal and elastic properties of the medium 
become temperature dependent, the formulas used before are 
no longer applicable. And finally, for significant values 
of efficiency, an appreciable fraction of the input energy 
is lost by radiation and so the temperature distribution 
assumed cannot be used. 

7. Thermoacoustic Arrays 

The way the conversion mechanism has been treated 
here through a simplified thermodynamic model emphasizing 
the process of conversion of laser light to acoustic energy 
in liquids, may be considered an engineering approach of 
great practical interest. 

This phenomenon can be analyzed in a more general 
sense to include the treatment of a wider range of specific 
cases with different types of media and electromagnetic 
sources. 

R.M. White [Ref. 10] in a theoretical and experimental 
work analyzed the conversion process for several temperature 
distributions due to different input energy fluxes. 

The generation of low frequency sound from the 
absorption of light has become a subject of growing importance 
with potential significance in underwater acoustics. 
Westervelt and Larson [Ref. 11] in a recent theoretical 
formulation, predicted the development of highly directive 
sound beams in water by means of thermalization of modulated 


laser light. The generation of this highly directional 
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acoustic wave results from an extension of the concept of 
a Simple thermoacoustic source to that of a thermoacoustic 
array, that is, a series of thermoacoustic sources where 
pressure fields combine to give a directional effect. 

A theoretical summary of the mechanism leading to 
the preduction of such a wave, which propagates perpendicular 
to the axis of the laser beam, iS given next. 

Assume that a laser radiates an amplitude modulated 
beam of cross sectional area So propagating in the +z 


direction as shown in Fig. 2. 





Fig. 2. Geometry of Laser Excitation 


Due to the absorption of energy in the medium, each 
elemental volume experiences a rise in temperature which, 
in turn, causes an oscillatory expansion in the thermal 
volume. If the beam illuminates the surface of a medium of 


low thermal conductivity, the solution of the wave equation 


a5 





for the pressure distribution resulting from the exponential 


absorption of energy is 


Faoarr | jw/v|r-5"| 
= OO: = you =O 2 e Se 
p the = fe a. ae ar 
|r-r' | 


where w is the laser angular modulation frequency, a the 
logarithmic coefficient of thermal expansion, C the specific 
heat per unit area, ry the output intensity of the laser, 
a the attenuation coefficient for the laser beam and v the 
sound velocity. This result was obtained for a source 
function representing the energy per unit volume per unit 
time absorbed, given by I vexp(-az-jut). 

Looking for a far field approximation, r has to obey 


Ene Condition 
ea (85) 


So, Simplifying and performing the integraion on eq. (84), 
the acoustic pressure field radiated by this mechanism 
turns out to be 


jawaP eJkr-jot 


Pp =“ ~ “Tire ° uw¥ykeos 6 ve) 


where Py is the laser output power. Averaging eq. (86) 


over one period, the resulting acoustic intensity would be 


given by 


56 





nwa Q) 2 1 
ac + roeoeae 


ao 


a ee ee f 
20C (Tare ey) 


In this equation ep represents the density of the medium. 
The maximum value of eq. (87) occurs at 6 = 71/2 
and the angle at which the intensity 1s reduced to one half 


its maximum is 


8 mo 


172 (88) 


bola 
I+ 
wie 


The total acoustic power generated when the limit condition 


k/a >> 1 applies, is given by 


2. 2 
awa EG 


at, ae (89) 
lL60C 


The radiation pattern given by eq. (87) 1s symmetric about 
the Z axis and has a sharp maximum at 6 = 7/2 in the limit 
Gees Ks 

Equations (88) and (89) show that when a increases 
giving rise to more intense thermal expansions, larger 
acoustic power outputs are produced, but when a decreases, 
the length of the array increases and so the directivity 
1s improved. 

On the other hand, for a given value of a, increasing 
the modulation frequency w, the beam may be made as narrow 
as desired. An increase in conversion efficiency is also 


achieved as w is increased. 


Sa), 





The array acts then like a broadside radiator 
where directivity is proportional to the attenuation 
coefficient. 

All these theoretical results were verified and 
validated by experiments performed by T.G. Muir, C.R. Cul- 
bertson and J.R. Clynch [Ref. 12]. These experiments showed 
further that for applications of laser excited arrays to 
ocean acoustics, high-power lasers inthe green-blue region 
are required in order to obtain extremely narrow acoustic 


beams. 


C. STIMULATED BRILLOUIN SCATTERING 
ia iloriltnocaleering 

Sound waves in a medium may also be created from 
the interference of two intersecting laser beams. This 
effect has been verified experimentally [Refs. 13,14] and 
its detailed theoretical formulation may be found in the 
literature [Refs. 15,16]. 

A scattering phenomenen can be defined for the 
purposes of the present discussion as a process whereby 
a physical system extracts energy from a beam of light and 
re-emits this energy at the same or at a different wave- 
length, in the same or different directions from the incident 
beam. 

Scattering may occur in a number of different 
Situations but among all the forms it can take, Brillouin 
scattering is the oneto consider as far as conversion of 


electromagnetic to acoustic energy is concerned. This type 
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of scattering is an interaction process between acoustic 
waves at microwave frequencies and electromagnetic waves 
at optical frequencies. 

PeiacsseootMeaceCenteaurnateal Clectromagnetic field 
has both a particle and a wave nature. Such a statement 
leads to the definition of photon as a particle charac- 
terized by zero mass, energy hy and wave vector k. Simi- 
larly, the treatment of lattice modes of vibration in 
material media, leads to the phonon concept, a particle 
having zero mass, energy hv and momentum hg Wileine vets 
the phonon frequency, g the phonon wave vector, h = h/2n 
and h the Planck's constant. 

The acoustic waves involved in the production of 
Brillouin scattering, usually range from ito LO Nome Hz 
and can be externally introduced into a material medium, 
can be already present in the material in the form of 
thermally activated phonons or can be stimulated in the 
medium by the light beam itself. 

Zee aggmnet lection 

Consider that a sound wave at angular frequency 
W is traveling upward in a cubic block of a material 
medium, as shown in Fig. 3 in a cross section perspective. 

When a plane wave of light of angular frequency 
Wo hits the left face of the block in a direction such 
that its wavefronts are perpendicular to the acoustic 
wavefronts, an interaction process will take place within 


the material. 


Do 
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FIG. 3 — The Debye-Sears Effect 


The results of this scattering process, called 
Debye-Sears effect, is similar to what happens when a 
carrier frequency We 1s phase modulated with a single tone 
Wa! that is, light at angular sideband frequencies Wor 
Wot Wot2u etc., will emerge from the right face of the 
block. If X and A represent, respectively, the light and 
sound wavelengths, then the wavefronts associated with the 
first upper and lower sidebands WotW ns will be tilted by 
an angle 9 ~ +i\A/A with respect to the undeflected carrier 


eaecy ate Wes Similarly, for the nth order sidebands, the 


deflection angle will be + nd/A. 
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The output intensity of the various sidebands is 
the same not only because their amplitudes have to obey 
Bessel functions behaviour, but also because a destructive 
interference occurs in the medium. For a fixed wa the 
extent in which each sideband is 'destroyed' depends 
upon its frequency as well as upon the physical dimensions 
of the material. So, this type of scattering will diffract 
just a small amount of the total input light. 

Suppose now that the input ray is rotated to make 
an angle a ~ A/2A with the acoustic wavefronts. The process 
developed is equivalent to that occuring when the incident 
ray was parallel to the acoustic wavefronts. 

As the first upper sideband still makes an angle 
Wes with the undeflected output light, the interaction pro- 
cess for this sideband may be though of as a simple reflection 
of light from the acoustic wavefronts. This effect is 


depicted in Fig. 4. 





(a) (b) 
FIG. 4 — The Bragg Effect Geometry 
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Moreover, this sideband is the only one suffering 
a constructive interference in the medium and therefore 
it 1s the only one with significant intensity. 

The angle of incidence a such that q = \/2A, or 


more precisely 


1 


F (90) 


a = sin 
1s called Bragg angle. Clearly, a Bragg angle can be 
defined for each output frequency. 

Due to this mechanism, the Bragg effect allows the 
selection of a particular mode of output light by a 
Simple control of the direction of the incident light 
that is, by forcing light and sound to meet at the appro- 
priate angle. 

It can be shown that at exactly the Bragg angle, 
the scattering attains its maximum. 

3. Production of Sound Waves 
When two light beams with frequencies Wy and w 


2 


-> ~ : 3 
and wave vectors Ki and Ko intersect in amedium, an 


MicucedmelectroStrIctive PmReSsure wave Of Lrequency 
|W Ws | and wave vector eke may be preduced- If the angle 
of intersection of the light beams is such that the phase 


Velocity Wy -Wo/ |K,-K Sree temGecslnemwave 1S Ene) Same 


2 | 
as the velocity of the sound in the medium, then a traveling 
wave will be produced. It is verified that the angle 


Satisfying this condition is precisely twice the Bragg 
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migue.) §£t snould be noted that this effect can only occur 
for an input beam energy above a well defined threshold. 

This mechanism of interaction of intense light 
beams leading to the production of coherent sound is called 
stimulated -Brillouin scattering. 

Consider now the situation shown in Fig. 4 where 
light and sound intersect at the Bragg angle. The inci- 
dent and scattered light wavefronts intersect at twice the 
Bragg angle and so they should form an interaction pattern 
qeving rise to stimulated Brillouin scattering. Thus, 
additional sound at Wa is launched from the interaction 
region. 

fie cumns Out that in the casevof Fig..4 (a), the 
produced (additional) and original sound waves are in 
opposite phases causing a progressive decrease of the 
sound amplitude (deamplification). On the other hand, if 
the situation shown in Fig. 4(b) occurs, the produced and 
Original sound waves are in phase and amplification of 
sound is then verified. 

The stimulated Brillouin scattering mechanism may 
be reduced to the following picture: if the input energy 
is high enough, instability conditions are created within 
the medium and, even when no input sound is present, pro- 
duction and buildup of sound at w, = Jw,-w,| starting from 
thermally activated phonons, will be verified. Or, in other 
words, the elctrostrictive mixing of two light beams gives 


rise to amplification of thermal sound in the medium. 
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Brillouin scattering is a classic example of 
parametric interactions in traveling wave systems where 
the power transfer between the waves may be calculated on 
the basis of either the Manley-Rowe relations or simple 
photon-phonon energy conservation arguments. It is impor- 
tant to realize that the case of stimulated Brillouin 
scattering gives an example where the principle of para- 
metric oscillations can perform a useful function, namely, 
the output energy at the lower frequency 1s in the form 
of sound as contrasted to the electromagnetic energy of 
the pump. 

It has been proved experimentally [Ref. 13] that 
the observed sound cannot be the result of thermal shock 
in the medium. In fact, removing one of the light beams 
without affecting the total input light power, it can be 
observed that sound is no longer produced. That is, the 
Signal observed at the beat frequency of the two light 


beams is indeed optically induced sound. 


64 





Preece ERI MENTAL PROCEDURE 


Pea Pe ARATUS DESCRIPTION AND RESULTS 
1. Radiation Pressure Experiment 

In the first experiment performed it was intended 
to observe the forces of radiation pressure developed in 
fresh water by a plane electromagnetic wave. 

Due to the mismatch between air and water, part 
of the energy is reflected back at the interface and part 
of it 1s transmitted through the water which is assumed 
here to be devoid of mechanical elasticity and of free 
charges. 

To perform the experiment, a portable radar model 
AN/PPS-6(XE-6) was used aS a source of electromagnetic 


field. Some characteristics of interest of this radar 


eeiecis 
BuUlcouscpeLttion rate —-----=----—--- a2 
ETS SY SAIN OY a a a ag oe aa Seo eGhZ 
Pete ( —————————= SSS SoS 02] sec 
ie a Ou Cage ne a a 100 W 
Veta Tmo OW Cig a= —————$— = Sa 47 mW 
power supply requirements ----------- ovat 2.58 


The power supply chosen was a POWER DESIGNS INC., model 
3240. Appropriate pieces of coaxial cable and rectangular 
waveguide both assumed lossless, were used to connect the 


radar output to a horn antenna, a SPERRY MICROLINE model 
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ae Lee AccOraGang tO 1lts calibration curves an absolute 
gain of 15.8 db can be expected at 9.3 GHz. The antenna 
was Oriented to face the surface of the water with the 
axils nearly perpendicular to it so that a transverse uni- 
form plane wave could be assumed at the illuminated surface. 

In order to detect the pressure waves produced in 
the water, an hydrophone type LC-10 was used. This is a 
small device about 5.4 cm long having a circular cross 
section with diameter about 0.96 cm. The useful steady 
state bandwidth of this acoustic detector extends from 
0.1 Hz to 120 KHz. The typical free-field voltage sensi- 
tivity within this range lies between -105 and -110 db 
relative to 1 V/ybar. The actual hydrophone used was the 
LC10 serial #2122 exhibiting a flat free field frequency 
response between 1 KHz and 9 KHz at -110 db. 

The free field voltage sensitivity, or microphone 


sensitivity level ML, is defined by the relation 


where V is the amplitude of the voltage output, P is the 
pressure amplitude of the acoustical field and 1V/ubar is 
the reference sensitivity. 

inewaelo typical directivity pattern in the plane 
perpendicular to the axis is almost circular. During the 
measurement, this device was held in the water in a posi- 


tion parallel to the surface and in the same vertical of 
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the transmitting horn antenna. The hydrophone depth could 
be varied since it was bonded to a horizontal plastic rod 
45 cm long which in turn was attached to a vertical 
metallic rod; this assembly could be moved vertically in 
front of a tape-measure stuck on a fixed wall. 

Through a 25 ft. coaxial cable the hydrophone 
Output signal was fed to a HEWLETT PACKARD voltage ampli- 
fier model 465A capable of either 20 or 40 db of amplification. 

Two different schemes were tried to measure the 
forces of radiation pressure. In the first, a cubic wood 
box about 6 ft. high was placed over one of the acoustic 
tanks in Spanagel Hall Room 025. The box was covered 
inside with a good absorber of electromagnetic energy; 
the transmitting equipment was placed on its top side 
where a small hole was made to let the horn antenna go 
inside. The hydrophone assembly was stuck to one of the 
vertical wood walls. The distance between the horn and the 
surface of the water was fixed at about 1 Mm. 

The output of the voltage amplifier was connected 
to the input of a TEKTRONIX oscilloscope type 422 either 
directly or through a very narrow band-pass filter. MThis 
filtering action was achieved through the operation as a 
narrow band-pass filter of a wave analyzer HEWLETT PACKARD 
model 302A. 

Due to the considerable environmental random noise 
amplitude and to the extreme weakness of the signal pro- 


duced by the energy conversion mechanism, no effect could 
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be observed in the oscilloscope. No signal component at 
about 2 KHz appeared when the output signal was passed 
through the band-pass filter, neither any signal resembling 
the envelope of the transmitted waveform could be seen 

over the noise when the output of the voltage amplifier 

was fed directly to the oscilloscope. The magnitude of 
‘the signal to be expected is computed later in the theor- 
etical computations.section. 

In the second scheme used it was intended to take 
advantage of a crosscorrelation or autocorrelation process. 

Instead of making the detection in the big acoustic 
tank mentioned before, a small portable water container 
was used. It was covered inside with styrofoam to minimize 
acoustic reflections; its internal dimensions were approxi- 
mately 26 x30 x67 cm. This time, the horn antenna was 
held at about 10 cm above the surface of the water. 

A tunable envelope detector was installed in the 
wave-guide in order to get the envelope of the transmitted 
waveform to be crosscorrelated with the received signal via 
hydrophone. A correlation and probability analyser SAICOR 
model SAI~43A with capability of performing either auto or 
crosscorrelation seemed to be a suitable equipment to use 
in front of the amplified signal. The output of this 
analyzer may be displayed in any ordinary oscilloscope. 

When the autocorrelation of the envelope of the 
transmitted Ta etocn was tried, no result was obtained 


since the detected pulse was so narrow that even in the 


68 








fastest sampling mode of the analyser no autocorrelation 
could be achieved on this waveform. So, instead of trying 
to detect the presence of the signal in the noise using 
PeesscOrrelation it was decided to make simply the auto- 
correlation of the noiselike waveform at the output of the 
voltage amplifier. The use of two voltage amplifiers in 
series became necessary in order to overcome a well defined 
input threshold level since only after this signal level is 
achieved the equipment can perform properly. 

A non-expected waveform was displayed in the oscillo- 
scope after the autocorrelation operation. Rather than 
being associated with the acoustic waveform expected, it 
was Simply the result of the autocorrelation of a 60 Hz 
power line interference corrupted with random noise. It 
could be observed that the same autocorrelation waveform 
was exactly obtained whether the radar was transmitting or 
not. 

As a conclusion of these experiments and according 
to the theoretical calculations, such a low input electro- 
magnetic energy is not enough to give detectable results 
after suffering this extremely weak conversion process. 

2. Transient Surface Heating Experiment 

In the second group of experiments performed the 
fresh water was assumed an elastic medium. If results were 
eventually obtained by using a high power pulsed laser 


they could be attributed mainly to a thermally shocked 
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interface rather than to the simultaneous but much weaker 
radiation pressure effects. 

Two types of lasers were used in this part of the 
experimental work as sources of electromagnetic energy. 

The first one was the Nitrogen pulsed laser existing 
in Bullard Hall Room 217. This device is capable of 
delivering in normal operating conditions as much as 80 
kw peak power. Its pulse repetition rate can be varied 
but 20 pulses per second was the one used. The pulse 
width is about 10 usec and the wavelength is 0.3371 u. 

The detection apparatus was again the hydrophone, 
voltage amplifier and oscilloscope. 

Unfortunately the output light of this source has 
not the shape of a narrow beam; instead it diverges too 
much and in order to focus as much light as possible several 
lenses had to be used. A mirror was also used to deflect 
to a vertical orientation the horizontal output beam. As 
a result of all these effects and characteristics the input 
energy at the surface of the water was not enough to generate 
a detectable signal at the oscilloscope display. Due to 
improper shielding of the laser, electromagnetic interfer- 
ence at the pulse repetition frequency was picked up also 
random noise and power line interference were observed. 

No conversion effect was then detected through 
this scheme. 

The second laser used was a single pulse laser 


existing in the Laser Laboratory, Building 223. This is a 
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high energy Neodymium Glass Rod laser operating at the 
wavelength 1.06 u and its output light is highly direc- 
tional. The beam was deflected to the vertical orientation 
uSing a low absorption prism. It is possible to vary its 
output energy and values up to 12 Joules can be safely 
obtained. However with the control knobs set at a fixed 
position it is observed that the output energy do not stay 
constant from shot to shot. The extent of this variation 
can be appreciated next in Fig. 6. 

Through the use of a beam splitter a portion of the 
output beam was diverged toward a photo detector whose 
output voltage signal was used to measure the energy of 
each output pulse. Another photo detector excited by the 
visible light of the laser flash tube was used for scope 
triggering purposes. This triggering action took place 
about 0.1 msec earlier than the beginning of the laser 
pulse. So, there is a time delay of about 0.1 msec between 
the starting of the sweep and the beginning of the laser 
pulse that should be considered when interpreting Fig. 6. 
Two TEKTRONIX oscilloscopes were used: a memo-scope type 
564B to measure the energy of each pulse and a 545A with 
Single sweep capability to observe and photograph the de- 
tected signal with a polaroid TEKTRONIX oscilloscope camera 
system. A sketch of the experimental setup is shown in 


Be Gis. 
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FIG. 5 — Experimental Setup 


The experiment was performed with sea water to 
decrease the penetration depth of the laser beam. In fact 
at pes x 10 Hz, which is the laser frequency, the pene- 
tration depth for sea water with productivity o = 4 mhos/m 
1s about 1.19 cm while for fresh water it is orders of 
magnitude greater. Thus, the use of sea water would guaran-~ 
tee a total absorption of the electromagnetic energy within 
the limited system used. 

The hydrophone was not exactly in the same vertical 
of the laser beam to prevent the possibility of damage 
caused by the intense radiation. 

To make sure that the signal obtained after each 
shot was not the result of electromagnetic interference, 


a flat piece of wood was placed over the container to prevent 


the radiation from reaching the surface of the water. 
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FIG. 6 — The output energy, water surface condition and 
hydrophone depth are shown for each shot. 
Horizontal scale is 0.1 msec/div for all pictures. 
Vertical scale is 5v/div except for picture (h) 
where it is 2v/div. 
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Picture 6(a) was taken under these conditions and as it 

can be seen no signal was recorded; this shows clearly 

that the signals obtained before came really from a surface 
heating conversion process. 

Fig. 6 presents the results of a series of shots 
for-various surface conditions and hydrophone depths. The 
voltage amplifier was kept at 20 db for all of them. As 
mentioned before although the laser controls were ata 
fixed position the output energy was not constant and 
values ranging from 5.7 to 7.1 Joules were obtained. 

From pictures 6(b) to (f) it can be seen that the 
time delay between the starting of the laser pulse and the 
beginning of the detected signal is about the same whether 
the hydrophone was at 5 or 10 cm depth and much bigger 
than the time that sound in water takes to travel 5 or 10 
cm. Obviously this is due to the facts that heating is 
not instantaneous and that the energy is not only abosrbed 
at the surface that causes the thermal effect and the sound 
production to take place but over a considerable region, the 
extent of which is determined by the penetration depth. 

Pictures (g) and (h) were taken with the hydrophone 
much farther from the beam vertical than in all the others 
and at 17 cm from the small surface area illuminated. Here 
the effect of a greater distance from the shocked region 
is apparent since the time delay is now greater. 

Pictures (d) and (f) were taken under constrained 


surface conditions. As predicted by the theory the signal 
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obtained had a higher amplitude; it also happened however 

that the energy in both cases was higher than in the free 

surface case and so it is not clear in which extent the 

increase in amplitude was determined by the surface constraint. 

It should be noted that a circular glass plate of about 

80 om? was used touching the water to constrain its surface 

which means that just a small percentage of it was constrained. 
Clear is the fact that in general the shape of the 

waveforms obtained is not close to the theoretical predictions 

made before. It is not difficult to accept such a different 

behavior since the assumed theoretical conditions were not 

fulfilled in this experiment, namely the observation point 

was not in the far field and the heating was not at or near 

the surface. Moreover, it is known that the envelope of 

the output pulse of this laser system is not a neat rectangle 

but exhibits instead two peaks and has a shape and duration 

varying from pulse to pulse. Also the relatively long plastic 

rod supporting the hydrophone may easily get into vibration 

when subjected to the acoustic pressure and acoustical 


reflections from the walls are likely to occur. 


TO appreciate the pressure field peak amplitude 
consider for example fig. 6(c) where a peak voltage amplitude 
of about 10 V was obtained. Taking into account the 20 db 
of signal amplification and using eq. (91) it can be 
computed that the pressure amplitude of the produced 


acoustical field is 3.16 bars or 3.12 atmospheres. 
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Pee LHEORETICAL COMPUTATIONS 

Theoretical calculations related to the radiation 
process are now made to give a further insight on the 
problem and to allow a better understanding of the diffi- 
culty in getting experimental results at low input power 
levels. 

Assume that the medium under consideration is sea water 
with conductivity o = 4 mhos/m, relative permittivity eee 80 
and relative permeability te le 

The reflection coefficient R for a plane electromagnetic 


wave incident normally on the surface of the water is given 


by 


elo 
Bi, ont ny 


where E- and E; are respectively the reflected and incident 


E field strengths, ie: = Mh Ve = S779. 1S the Intrinsic 


impedance of free space and 


= =) oe ee 
Cunds JWE LEY 


is the intrinsic impedance of the water. When the frequency 


of the radiation is £ = 9.3 GHZ then 


2 


eS mcit tooo LO mrad /7 Sec? 


n = 42.02 Zo Os; 
R = 0.799 Oto S : 


76 





This means that only 20.1% of the E field will penetrate 
the water. 
If the incident electromagnetic energy is from a radar 


with output peak power Py = 100 W and if a horn antenna 


with gain c= 15.8 db = 38.02 is at r = 1m above the 


surface of the water, then the power density at the interface 





is 
2 
P.G 5. 
ea Ss w/t ee 
a 2 n 
4nx Oo 
215 Es =n, * 302.55 or E, = 337.73 V/m. Then 


Oy See. 1 29 — 6 oo 7m 


will penetrate the water. 


The attenuation constant of the water is given by 


WZ 
ie Cha = 
a = ol 5 { at re ae = 84.142 Np/m 


The penetration depth 6 = 1l/fa is then equal to 1.188 om. 
So, as far as absorption of energy 1S concerned, the medium 
may practically be considered as extending from z = 0 to 

Z =o, Thus the components of the radiation pressure force 
are in this case the force given by eq. (46) due to energy 
dissipation of 20.1% of the incident energy and the force 
given by eq. (48) due tothe reflection at the interface of 


79.9% of that energy. 
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From eq. (46) the average force due to dissipation of 


the transmitted energy is 





ae 2 ) cos 6 
O q 
maere |yn| = 42.02 2, 6 = 2276 and Eg ~ 67-88 V/m; v is given 


by eq. (24) and = Wye 


Computing the values of v and ve Teeeurno GUL that 
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11.23 x 1077 . 
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The average force is then 


8 Z 


163.195 x 10 ~*~ N/m 


rr 
I 


U63,095 x 107? dyne/em? 
7 


163.195 x 10. 


ubars. 


With this value for the average pressure field and knowing 
that the hydrophone sensitivity level is ML = -110 db 
relative to 1V/ubar, it is possible using eq. (91) to compute 


what would be the hydrophone output voltage 


Sel 0 


20 log V/Fy 


log V/Fy -5.5 
V/F, = 3.16 x 107° 


aCe creonlon 


Volts. 
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From eq. (48) the average force due to energy reflection 


at the interface is given by 


where the factor ve UL 1s now 1 since incidence and 
reflection of energy are taking place in free space, 


c= 3.x Loe m/sec and 
2 
s = = = a = 96.577 W/m 


SO 


8 2 


64.385 x 10 ~ N/m 


Bice een eee 


Following the same procedure as before the output voltage 


would be 


0.203 x 10°79 volts. 


~*~ 
t~ 
© 
t 
Il 


The values just obtained for the output voltages due 
to both force components are so small that they cannot be 
easily seen in ordinary oscilloscopes. More sophisticated 
detection schemes could lead to positive results but as 
mentioned before even through a correlation process, it was 


not possible to detect the presence of the signal. 


og 





Electromagnetic energy has to be increased orders of 
magnitude if radiation pressure forces are to be easily 


displayed by detection equipment. 
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ive CONCLUSTEONS 


According to theoretical predictions and experimental 
results, in order to get significant acoustic power, 
extremely high electromagnetic energy has to be used. The 
reason for this is the low efficiency that characterizes 
this type of energy conversion. 

In practice, when the conversion process takes place, 
both radiation pressure and transient surface heating forces 
will be developed in the medium. 

Thermal forces are by far more important since they 
are orders of magnitude larger than radiation pressure 
forces. The present status of laser technology has spurred 
interest in the possible applications of this conversion 
process. In the future, it is likely that there will be 
ocean applications of highly directive thermoacoustic 
arrays excited by modulated lasers. 

Radiation pressure forces, although much smaller in 
magnitude and requiring careful design of detection schemes, 
higher power sources and highly sensitive devices, have been 
proposed to be used in several applications. 

A. Ashkin [Ref. 17] observed that micron-sized particles 
in liquids and gas may be accelerated and trapped in a laser 
beam. The experiments were performed taking special care 
to avoid as much as possible thermal effects from the laser 


light so that the observed effects were due primarily to 
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radiation pressure forces. Chemical and laser pumping 
applications were suggseted on the basis of this effect. 

H.E.M. Barlow points out in his paper [Ref. 4] the 
possibilities of the application of radiation pressure 
forces in power measuring instruments and in orientation 
behavior in satellites. 

" Much work has to be done yet in the field of the 
conversion of electromagnetic to acoustic energy until all 
useful applications are thoroughly explored. 

As a final comment it should be noted that the Laser 
Laboratory at NPGS seems to have the appropriate material 
available to encourage further experimental work on the 


production of acoustic waves in media through transient 


surface heating. 
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APPENDIX A 


DETAILED THEORETICAL CALCULATIONS 


Penetration depth (6) 


for sea water 


o = 4 mhos/m 
e = 80 
= Sy: 
ce = Ef r =e. 32 lO F/m 
Ee = 1 


w = 2rf= 2 x 9. 

meee an 
= 8.9 

noe. 

WE 58.43 x 107 


a = w = ( Ji=e 


58.43 x 10° 


58.43 x 10° 


84.142 Np/m 


020168" m 


Q(t 


For the radar frequency f 


= 9.3 GHZ 
9 ¥ 
3.x 10 >= 56545 210° rad/sec 
x oma x 8.854 x 10 14 x 80 
peo ee 
4 ee Le 
a2 eee 


rel. 1) 


/ 4.4505 x 10 °° x 0.0046596 
/ 0.020738 x 10 7° 


2S 6 cem 
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= 0.0966484 


Ea 10 ny AO emo 0oosai 9 = 91) 





For the laser frequency 


8 
fees A, = 2.83 x 1014 Hz 
ieO6 x 10 
_ _ 14 
jee weonte —wie. 78 x LO rad/sec 
O 4 4 


in 12598902206 


Ceo 


17.78 x 1074 / 4.4505 x 10 +© x» 5 x 107° 


Qa = 
= 83.873 Np/m 
ee = =eOnOL oem = latoomen 


Reflection coefficient (R) for sea water for the radar 
frequency 


= 
li 


i 


a LS 


cr 


58.43 x 10° rad/sec 
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1 = | —JeY _ | ji58.43 x 107 x 47 x 1077 
+ 7 ea. Qo ea 
ce 4 + 458.43 x 10° x 708.32 x 10 2° 


_ | TBST ees sree ee 
{/-4 + 441.3071 11579903 87248 
: ° 
AVEO Aue = 47974 ae 028 5 


ae vot eno coe — 30a = 3 e026 2.0235 
aneomeet 72.0235 + 377 4718-97457 32.0225 


9eI0°3461 


D°D7G7 0.799 |-0°623 


335.03 
418.979e2 





Radiation pressure force Fy due to energy dissipation 


= 
—- Vv xO 
Ss re x Dia o°° 8 
O 


— a a a ee 
Vo we[Yl + (a/uc)* + 1] 


2 


8.901 x 102° ,/ 1 + (0.0966484)° + WJ 


2 


VT = 


MIE 


= 0.3348 x Gi m/sec 





KG 


8.901 x 10 -° x 2.00466 


v2 = = ——— = 0.11235 x Oia 
- 8.901. x 10 


8 2 
{= es + 67.887 x 0.9988 = 163.195 x 10 
OL235 x 106 2. R24 22602 


F 
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